A Generalization of Grover's Algorithm 



Luigi Accardi 
Ruben Sabbadini 

Centro Vito Volterra 
Universita degli Studi di Roma "Tor Vergata" 
Via Orazio Raimondo, 00173 Roma, Italia 
accardi® volterra. mat .uniroma2 . it , |http: / /volterra. mat. uniroma2.it 



ABSTRACT 

We investigate the necessary and sufficient conditions in order that a uni- 
tary operator can amplify a pre-assigned component relative to a particular 
basis of a generic vector at the expence of the other components. This leads 
to a general method which allows, given a vector and one of its components 
we want to amplify, to choose the optimal unitary operator which realizes 
that goal. Grover's quantum algorithm is shown to be a particular case of 
our general method. 

However the general structure of the unitary we find is remarkably sim- 
ilar to that of Grover's one: a sign flip of one component combined with a 
reflection with respect to a vector. In Grover's case this vector is fixed; in 
our case it depends on a parameter and this allows optimization. 



Unitary operators which increase the prob- 
ability of the |0 > component of a pre-assigned 
vector 



Let \i > (i — 0, . . . , N— 1) be an orthonormal basis of R . The mathematical 
core of Grover's algorithm is the construction of a unitary operator U which 
increases the probability of one of the components of a given unit vector, 
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in the given basis, at the expence of the remaining ones. The necessity of 
such an amplification of probabilities arises in several problems of quantum 
computation. For example in the Ohya-Masuda [4] quantum SAT algorithm 
such a problem arises. In a recent interesting paper Ohya and Volovich 
have proposed a new method of amplification, based on non linear chaotic 
dynamics [14]. In the present paper we begin to study the following problem: 
is it possible to extend Grover's algorithm so that it becomes applicable to 
a more general class of initial vectors, for example those wich arise in the 
Ohya-Masuda algorithm? A preliminary step to solve this problem is to 
determine the most general unitary operator which performs the same task 
of Grover's operator. This is done in Theorem (1.1) below. The result is 
rather surprising: we find that, up to the choice of four ±1 (phases), there 
exists exactly one class of such unitary operators, labeled by an arbitrary 
parameter in the interval [0, 1]. Moreover these unitaries can be written in 
a form similar to Grover's one, i.e. a reflection with respect to a given unit 
vector possibly preceeded by a sign flip of one component combined with, 
where the unit vector in question depends on this parameter in [0, 1]. The 
free parameter in our problem allows to solve a new problem, which could 
not be formulated within the framework of Grover's explicit construction, 
namely the optimization problem with respect to the given parameter. We 
prove that, even in the case of Grover's original algorithm, this additional 
freedom allows to speed up considerably the amplification procedure. In a 
forthcoming paper [15] we plan to apply the present method to the Ohya- 
Masuda algorithm. Since an operator U is unitary if and only if it leaves 
unaltered the scalar products of vectors with real components in a given basis, 
we shall restrict our attention to unitary operators with real coefficients in a 
given basis (as the Grover's ones). This restrictes the problem to R N . 

THEOREM 1 Given the linear functionals: 



f] : a = (di) e R N h-> 77(a) = Vi a i 



(1) 



i=0 



N-l 

c: a= (di) G R N 1— > c(a) = Y 7i°i 

i=0 



(2) 
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with 7j and rji real and e±, £2 £ {±1}; necessary and sufficient condition for 
the operator U, defined by: 

U^2cii\i >= Ei(a + r](a))\0 > +e 2 ^2 ( a * + c ( a )) N > ( 3 ) 

to be unitary is that there exist a real number (3q such that: 

IA>I < 1 (4) 



7, =- 1 ^r ^° (6) 

r] =-l + e A /3 (7) 
?7i = ^37o (8) 
where £3, £4, £5 are arbitrarily chosen in the set {±1}- 

PROOF In finite dimension unitarity is equivalent to isometry. Therefore U 
is unitary if and only if, for every \a >= X^Lq 1 a *K > the following isometricy 
condition is satisfied: 

= (a + r]) 2 + Y,(ai + c) 2 = a 2 + if + 2a r] + J2 a l + ( N ~ l )<? + 2c E a i 

where we write 77, c for 77(a), c(a). This condition can be written in the form: 
77 2 + 2a T7+(iV-i)c 2 + 2c5> i = (9) 

With the notation: 

7 (a)= 7 :=(iV-l)c 2 + 2c5> i (10) 

Equation (®) is equivalent to: 

t; 2 + 2a 77 + 7 = (11) 
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and its possible solutions are: 

77(a) = 77 



-ClO + B4\Joq - 7(a) 



(12) 



Given (p~2|) the funtional 77(a) will be linear if and only if Va , 



a N : 



a l ~ 7(a) = E P: 



3^3 



(13) 



for some real numbers j3j indipendent of a. 

Since the functional c(a) is linear and given by (0), because of (|12|) and 
(|13"D, condition (|9|) becomes: 



-aS + (JV-l) 



2 ' \ 2 

E [,<>:, ■ -E7/'oE"< 



-al + 2 E Ijdj J2 a i + J2l( N ~ l)7t7j + A&] a * a j = 
or equivalently: 

(N - l) 7o 2 + $-l] + E t 2 ^ + (JV - l)7i7i + PiPj] Wj+ 

+ 2 E [7o + (iV - 1)707* + A) A] aoa* = 



(14) 



The identity ( |14[ ) holds Vao, . . . , a at, if and only if: 



(iV - l) 7o 2 + (3 2 - 1 = (15) 

2 7 . + (AT -l)jrf3 + PiPi = Q Vi.i^O i^j (16) 

2 7l + (N - 1) 7 2 + /3 2 = V7^0 (17) 

7o +(iV-l)7o7, + A)A = Vz^O (18) 

Equation (|U|) and the reality condition on rj imply that (|J) and ([5]) hold. 
From (^) we deduce that, for % 7^ 0: 

70 + A)A 



7i 



7o (iV-l) 



(19) 
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and, replacing this into (|T7|)> we find: 

2(7o + A) A) , (7o + A) A 



\2 



+ w ?/,r'~:( +a 2 = q 



7o (iV-l) 7o 2 (^-l) 



or: 



(iV-l)7o 2 + ^ A" = 7o 



which, because of (|15f) , is equivalent to: 



l-/?o 2 



A = £s7o = £3^5 \ N J ° (20) 

with £3 = ±1. Replacing (|20|) into (|l^) we arrive to (|B]) 

Replacing (0),. . ., (P^ ) into (]!]) and (0), we conclude that a necessary- 
condition for the linearity of U is that 77 and c must have the form: 



v^A 



77(a) = (-1 + £ 4 A))«o + £4^370 51 a ^ = ( _1 + £ 4/?o)ao + e A e 3 e 5 ^-j== ^ a fc 

fc^o V N — 1 fc^o 

(21) 

c(a) = 7 a - — — — 2^ a fc = 65^==^ - — — — ^ a k (22) 
iV ~ 1 mo viV - l iv - 1 M0 

Conversely, if conditions ([|), . . . , @ are satisfied, then also (0), which 
is equivalent to (||), is satisfied and therefore U is isometric, hence unitary. 
This can also be seen by a direct computation (see appendix A). 

REMARK Because of (|]) there exists a # G [0, 2n) such that A> has the 
form: 

A) = e 3 cos 9 (23) 

and therefore, from (^|): 

VN^llo = e^l-ffi = sm 9 (24) 

i.e. the parameters A> an d 7o live onto an ellipse in the (A), 7o)-plane. With 
these notations one has: 

SITl 9 

77(a) = (-1 + e 3 e 4 cos 9) a + e 3 e A y= — = ^ a k (25) 

V A - 1 fc ^ 
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sm 



c a 



a 



1 + cos 6 
N — 1 



k^O 



(26) 



REMARK The case 7 = leads to r] = or i] = — 2ao; in both cases we 
have: 

U }^aj\i >= ±£ia |0 > +52 XX a * + c ) I* > 

The operators C/(7 = 0(a)) are in this class, however they play a significant 
role in Grover's algorithm because they may be used to change the sign of a 
component leaving the others u inaltered {flip)- 

If we are interested in unitaries which modify the component of a, we 
must look for solutions with 7 7^ 0. 

COROLLARY 2 If in fl2I|) and we choose: 

6164 = £3 = £5 = 1 
e 2 = -l 

2 

l0= N 

then the corresponding operator U is Grover's unitary (see section 4). 

PROOF It is known that Grover's unitary is characterized by (see section 
4): 

N — 2 2 ^ 
«o >-> — ^r~ a o + T7 2^ a k ='■ £ i l a o + VW\ ( 27 ) 

a* -a* + [ -a + 51 a fc J =: £ 2 [oj + c(a)] (28) 
On the other hand, from equations (|25| ) and ( p£j) we have: 

£1 [oo + 77(a)] = £1^4 /?o«o + £370 a k ( 29 ) 
V Mo / 
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e 2 [Oj + c(a)] = e 2 \ai + 7 a - 1 ^^ 3 ^° X] a ^ 



with 70 given by (|5|). Comparing this with ( 
condition for equality is: 

R - N ~ 2 
£\£aPq — — — 



and 



Now let us choose E1E4 = 1 and /3 = ^j^- then: 



7o = £5* 



1 



(JV-2) 2 



\ N — 1 



£5 N 



(30) 

we see that the 



that leads to £5 = 1. Therefore, if e 2 = — 1, the coefficient of the third term 
in fl30T) becomes: 



-e 2 - 



1 + e 3 A) 1 + £3 



iV-1 



jV-2 
N 



N-l 



N + e 3 N - 2e, 
N(N — 1) 



that gives the correct parameter if and only if £3 = 1. 



2 Canonical form and reflections 



THEOREM 1 Any unitary operator U(/3o,e) with real coefficients in the 
basis (|z >) and satisfying the conditions of Theorem (1.1), can be written 
in the form: 



U(j3 Q ,e):=e 1 e 4l \0> |/3 < 0| + e 3 e 5 ^== £ < k\ J + 



k^O 

Moreover a unit vector u G R N such that: 

U{fa,e)=e 2 {\-2\u><u\) (32) 
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exists if and only if e is such that 

£2 = £i£i£3 

In this case \u > has the form: 



u > 



1 

V2 



-£5 



r 7— / l + e 3 A ) 



'1+eaA) 
iV — 1 



(33) 



(34) 



REMARK Notice that unitary operator fl32|) simply realizes the reflection 
of the \u >-component of any vector with respect to the \u >-axis. 
PROOF The identity follows immediately from (@), (^TJ) and (|22|). 

The operator U(fl ,e) of the equation ( |31~D can be rapresented in the 
following way: 



U{(3o,e) :=£ 2 1-<M0> 



+ e 2 J2 N > 



(e 2 l - £i£4A)) < 0| - ^i^^s 



+ 



-£5 -^7 5 < °l + — — r E I 



iV 



(35) 



Now an easy calculation shows that, given a vector \u > of the form (|35|). 
the right end side of (|3"2|) is equal to: 



e 2 l 



10 > 



+ e 2 S|i > 



£ 2 (1 - £3^0) < 0| - £ 2 £5 



1 £fn 



" 5 Va^T < 01 + iv 



1 + £ 3 /?o 



+ 



1 fc^o 



(36) 



For (3q 7^ 0, 1 (|36| ) and fl35|) are equal if and only if ([33]) holds, and the 
last operator is a projector if and only if: 



£2 — £\£a£^ 



(37) 



because the off-diagonal terms must be equal. From this the thesis follows 
observing that multiplying fl3"T|) for e 2 ^3 we obtain: 



£3 — £\£a£2 



(38) 
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COROLLARY 2 Grover's operator is the product of a operator of the form 
(32) with a flip, realized with a operator {7(7 = 0, Va), as in the Remark 
after previous Theorem 2.1. 

PROOF As it is implicit in its definition (see Section 4), Grover's operator 
is a flip followed by a reflection of the \v >-component with respect to the 
\v >-axis, where \v >:= N~ l l 2 \l, . . . , 1 >. 



REMARK Theorem (2.1) shows that Grover's unitary and the generalized 
ones presented in this paper are analogue, and the realizability of the former 
implies the realizability of the latter. 

REMARK If in (R5|) the identity fl3~T|) holds then, remembering (B3) and 
(|24|) , we can rewrite (|35|) inthe form: 



u > 



9 

— sin - 



cos I 



COS 



2' v/AT 



(39) 



which, up to a phase, is the most general form of a vector in R with N — 1 
components equal. 

REMARK A matrix rapresentation of the operator U(f3o,e) in the basis 
(|i >), with £ 2 = £15463 is: 



/ 



U(J3o,e) 



£2 



1 + £ 3 /3 - 1 
£5 





50 


v^v 


-1 




% 


VJV 


-1 








-1 



£5 



1 - 



^JV 11 ! 

l+£3/3p 

JV-1 

l+£ 3 /3 

JV-1 



l+£3/?0 

JV-1 



65 ViV-1 

l+£ 3 /3 
JV-1 

1 _ 1+£ 3 A) 
1 JV-1 



1+£3<3q 
JV-1 



fc 5 ^^zj 

1+£3/3q 
JV-1 

l+£3/3o 



JV-1 



I _ 1+E3A) 

x JV-1 / 



£ 2 1 - ^2 



-£5 



V -£5^3 



-£5 



1+£ 3 /3q 
JV-1 



l+£3/3o 

JV-1 



-£5 



VJV^T 

l+£3/3p 
JV-1 



l+£3/3o . 
JV-1 / 
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(i- cose -f^t 

sin 9 l+cos 9 



sin 9 



e 2 l - e 2 



sin 



e 2 \ - e 2 2 



sin 



__ sin 9 



JV-1 



1+COS I 

JV-1 



sin t^cost 



JV-1 



JV-1 

= e 2 (1 - 2\u >< u 
with \u > given by (p5|) or (|39|). 



x/JV- 

1+COS I 

N-l 



l+cos 9 . 



VN^l 



JV-1 



JV-1 



(40) 



3 Optimal Choice of the Parameters 



In this section we study the following generalization of Grover's problem: 
given a fixed vector \a >= J2i a i\i >> we l°°k f° r a unitary operator U = 
U(Po,e) of the form discussed in sections (1) and (2), which increases the 
probability of the 0-th component of \a >, i.e. such that: 

|oo| < \(U(Po,e)\a>) \ := \< 0\U(p ,e)\a >| (41) 

DEFINITION 1 A unitary operator U((3q, e) of the form discussed in sections 
(1) and (2) is an optimal amplifier for the 0-th component of \a > if condition 
(flTD is satisfied and: 

|(Z7(A),e)|o>) | < \(U(p ,£)\a>) \ (42) 
W(3 e [0, 1] ; := (ei, . . . , e 5 ) G {±1} 5 . If moreover: 

\{U(P ,e)\a>) \ = l (43) 
then we speek of an absolute optimal amplifier. 
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THEOREM 2 Given a unit vector of the form: 



\a G >:= a \0 > +b^2\i > (44) 

with ao 7^ and 

<2q + (AT — l)b 2 = 1; (45) 

an absolute optimal amplifier exists and is defined by: /5q = ±ao- 
PROOF From equations (gift, (H), (H) and © we have: 



U\a G >:= U I a |0 > +6 V \i > 



+£ 2 



(J a + e 3 e^(N-l)(l-Pi)b 
b + e 5 



10 > + 



l-/?o 2 



a - (1 + £ 3 /3 ) & 



/ / \ / -sen \ ^-^ 

£16453 (^cos 6>a + v — lsen 9 bj |0 > +£2 ^ ^/^y ~~ cos ^ / ,N 



> 



The amplitude of |0 > is extremal if: 






— — (cos 9ao + y/N — lsen 9 b) = —sen 9 a$ + VN — lcos 9 b = 
and this is satisfied by a such that: 



tg 9 = ^W^l — 

a 



(46) 



that gives: 



and 



£3/^0 = cos 9 



£6 



£5 vl — Po = sin 9 



1 + tg 2 9 

e 7 tg 9 
1+tg 2 9 



e^ao 



e 7 VN -lb 



where £ 6 , £ 7 G {±1} and we used (p>|). From ( f4~6"D we have £ 6 = £ 7 . 
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Therefore we obtain: 

a i — ► 6164638$ (sen 2 9 + cos 2 O^j = 81646385 

and _ _ _ _ 

/ sen Ocos 6 sen 6 cos 6\ 
1 > e 2 e 6 — / / = 

Thus the extremal amplitudes correspond to probability 1 and this completes 
the proof. 

REMARK The absolute optimality of the previous Theorem (3.2) refers to 
the case when all the components a,k(k 7^ 0) are equal. However, for a general 
vector, an optimal amplifier will not be absolutely optimal. This fact will be 
apparent from the following theorem. 

THEOREM 3 An optimal amplifier for a generic vector a of the form: 

N 

k>:=E a ili> ( 47 ) 
3=0 

with: 

5>; = i (48) 

exists if J2k^o a k 7^ and it is given by an operator (|31|) with the following 
choice: tg6 = 

PROOF From equations (ETh, (E3), pg) and <M) we have: 



N 



l-/?o 2 



U\a >:= Uj2 a j\J >= £ i £ 4 + £3^5 ^7= — = a k ] 1° > + 

3=0 V * ~ 



i 1 + 8 3 p \- 



\i >- 



j _ sen 6 m-^ \ 1 . 
£1^3 ( cos 6>a + 2^ a fc I 1° > + 

V M - 1 fc^O 
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^ / sen 9 1 + cos 9 ^ \ 

and the amplitude of |0 > is extremal for: 



% >= 



d ( sen 9 „ \ 

— \cos 9a + = = 2^ a k J = -sen 9a + 



cos 9 



then for a 9 such that: 
this gives: 

e 3 p = cos9 = — — 
V 



t g9 = 



E & ao\f N — 1 



l + tg*6 ./i + gy B f v /ag(iV-l) + (E M o 

v^T 7o = sen = = , £ ^^ ak 

^l + tg*9 v /ag(iV-l) + (E^o^) 2 



a k ) 



with e 6 , e 7 G ±1 and e e e 7 = 1, i.e. e 6 = £7- This gives: 

£i£ 4 £3£6 Oo(-^ - !) + (Efc^o «fc) : 
y/N=Ty/a 2 (N - 1) + (E^o «fe) 2 

e^e^Jc^N - 1) + (E^o Ofc) 2 

Finally: 



+ 



Efc^O a k) 

N - 1 



(49) 



^6 



ViV-l v /a2 (A r_ 1) + (EMo 



=o - 



^ _|_ g. apViV-l 

6 ^al(N-l)+(E^o a ^ 2 



N-l 



k^O 



e 2 \a i - 



Z^fc^O a k 

N- 1 



(50) 



and this completes the proof. 
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REMARK Obviously (g9]) and Qjjg) shall be as in Theorem (3.2) if the vector 
\a> is of the form ([44]). 

REMARK The action of the optimal amplifier found in Theorem (3.3) can 
be described in the following way: "For every \a >: 

1. we subtract from every a^, i ^ 0, the average of all the components 
different from the 0-th one: 



2. Then for the 0's component we have of course: 



Of) 



i — ( dj 



J2k^0 a k 

N — 1 



\ 



1-J24-(N-1) 



(Efc^O a k) (Si^O a i) (Sfc^O 

' (N-iy + 2 ~ 



a k 



N-l 



Y 0+ N-l 



(where in the last step we used PS])) as in the fl4"^)". 



4 Grover's algorithm 



Grover, in [1], considers the following problem (cf also [2]): 
PROBLEM: Given a (quantum) system with a state space 7i of dimen- 
sion N = 2 n . Let {0,1}^ = {So, S2, . . . , SV-i} = : S be the set of states 
represented as n q-bit string G Ti. Let be given a function: 

C: {0,1} 
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with the following property: there is only one state, say S v , such that C(S V ) = 
1, while C(S) = WS ^ S v . Construct a quantum computer algorithm 
which is able to find the unknown S v state with high (say > .5) probability. 

It is always possible to rename the states so that {Si, . . . , Sn} — {0, . . . , N- 
1} and S v = 0. In these notations let be given a vector: 



The first step in Grover's algorithm is to construct an operator Z that 
flips the 0-component. In our notations: 



Z : = 1-210 >< 01 



Grover then defines: 



\a >:= Z\a >= — a |0 > + ^<2j|z > 
and chooses the unit vector in formula (B2T) as follows: 



1 v-i, 1 h 

\v >:= — t= 2^ \k >= — ^=|1, . . . , 1 > 



N 



N 



(51) 



This gives: 



< v\a >= ^/=H < k \ ( — «o 1 > + ^2ai\i > J = ( -a + Y a k 



Then, denoting P := \v >< v\, Grover introduces the unitary operator 
D\a >:= —1 + 2P, whose action on \a > is given by: 



D\a >:= (-1 + 2P) \a >= —\a > +2 < v\a > \v >= —\a > + 



N 



"«0 + CL k 



a ° + ^ ak 



|o>+E 



\l > 



Then: 



N-2 2 _ 
a i * — — a + — 2^ a fc = £ i [«o + 77(a)] 



(52) 
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en h-> — a, + 



N 



-a + J2 a k) = £ 2 [oi + c(a)\ 



(53) 



If afc = Vfc, h ^ (the Grover's agorithm case) then: 



a 



N-2 2(N-1) 
-a H - a, ; 



JV 



2(N- 1) 



JV 



We can arrive to the same result working only via operator algebra. 
Grover's unitary is therefore: 

U G := DZ = -(1 - 2\v >< v\){l - 2|0 >< 0|) = 

4 



= -l + 2\v >< v\ + 2\0 >< 
and it acts on \a > in the following way: 

(-l + 2\v >< v\ +210 >< 01 - 



N 



v >< 



N 



\v >< 0\)\a >- 



N-l 



N-l 



"0 + E a J - M a °> I ai + AT E fl J ~ ^?°0 



j=0 



iV 



iV 



j=0 



JV 



i=l,...,iV-l' 



l-|)«o + ^E^ |«o-a, + ^E^ 

J V MO / i= i...JV-l' 



(54) 



Comparing (Fjl) with (f4D|)and ( 150 ) we conclude that, within the class 
U(f3o,e), Grover's unitary is not optimal, not only for a generic initial vector 
| a >, but even for the initial vector \ac > used by Grover himself and given 
by (fD. 
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5 A One Step Solution for Grover's Problem 



In the notations of section (4) let us write \j > for the state Sj so that, in 
particular, \S V >= \v >, and suppose moreover that j = 0,1, . . . , N — 1. 
In Grover's algorithm one uses the unitary operator: 

V\j>= (-l) c(s %> 

which is a self adjoint involution, i.e. V = V* and V 2 = 1 
We will consider the unitary operator: 

if j ^ or C(Sj) = 
V\j>={ |0> ifC(Sj) = l 




if j — an d C(S V ) = 1 

which is also an involution (V 2 = 1). 

Notice that V\j >= J2k u jk\k > and Ujk = Sjk if j ^ v, 0; u vo = 1, u v k = 
for k ^ v, u 0v = 1, u k = for A; ^ v. Thus V is a local operator in the sense 
of Grover [1], section (5). The action of V can be compactly described by: 



>= (1 - ^l,0(fi , o))*0,0(S,)b" > +£l,C(Sj)|0 > +(1 - 5l,C(S ))^',0+ 

+ ^E[l-(-l) C(5ft) ]l^> (55) 
z fc=i 

from which it is clear that the physical action of V is realized by parallel 
computation of the values of C, exactly as in Grover's algorithm. 
THEOREM In Theorem (3.2) let us choose the initial vector \og > in (fl5|) 
so that ao = b = -j= (i.e. we choose Grover's initial vector) and let us denote 
U the corresponding absolute uptimal unitary operator. Define: 

Uqpt :=V*UV 

then Uopt is an absolute optimal amplifier for the component > of \clq >. 
PROOF Clear from theorem (3.2) and the definition of U. 
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6 APPENDIX A: Direct proof that U verifies 
the Isometricity Condition 

Let us now verify that the set of conditions ([|), . . . , (||) are also sufficient 
conditions. To this goal we check if the isometricity condition ([ll|) is satisfied 
by the operator U, given by (0) if the parameters satisfy . . . , (|8|). Then, 
replacing (|2l|) and (H), which are equivalent to (|J), . . . , (|]), into (|T0|) and 
the (|12|), we find: 
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Therefore the isometricity condition ( |TT]) is equivalent to: 
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or equivalently: 
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which is an identity for any choice of £3, £4, £5 and this ends the proof. 
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